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Deterministic rate equations are widely used in the study of stochastic, interacting particles 
systems. This approach assumes that the inherent noise, associated with the discreteness of the 
elementary constituents, may be neglected when the number of particles N is large. Accordingly, it 
fails close to the extinction transition, when the amplitude of stochastic fluctuations is comparable 
with the size of the population. Here we present a general scaling theory of the transition regime 
for spatially extended systems. Two fundamental models for out-of-equilibrium phase transitions 
are considered: the Susceptible-Infected-Susceptible (SIS) that belongs to the directed percolation 
equivalence class, and the Susceptible-Infected-Recovered (SIR) model belonging to the dynamic 
percolation class. Implementing the Ginzburg criteria we show that the width of the fluctuation- 
dominated region scales like N~ K , where N is the number of individuals per site and « = 2/(d u — d), 
d u is the upper critical dimension. Other exponents that control the approach to the deterministic 
limit are shown to depend on k. The theory is extended to include the corrections to the front 
velocity above the transition. It is supported by the results of extensive numerical simulations for 
systems of various dimensionalities. 

PACS numbers: 87.23. Cc , 64.70. qj, 05.45.Xt, 05.45.-a 

INTRODUCTION 

The connection between a stochastic model of particle reactions (or equivalently, birth-death processes) and its 
associated deterministic rate equations is a topic of continuing interest. The common intuition is that the rate 
equations are not only qualitatively correct, but indeed provide, when the number of interacting particles is large, a 
quantitatively accurate approximation. This intuition is given concrete support by the Q expansion of van Kampen pQ . 
However, there are a number of situations in which this picture is too naive and needs to be refined. One by now classic 
example of this is the exponentially small rate of extinction for a system with an absorbing state [5] , which dominates 
the long-time dynamics, and is completely missed by the rate equations. Another example is the anomalously large 
corrections [31 14] to the front velocity in stochastic systems which exhibit propagation into an unstable state; e.g., 
systems whose rate equation is the Fisher-Kolmogorov equation. 

A system which captures features of both these examples is the spatially extended version [5] of the classic SIS 
(Susceptible-Infected-Susceptible) infection model of Weiss and Dishon [3] . In this model, contact (either on-site or 
nearest-neighbor) between an infected individual and a susceptible can, with some probability, convert the susceptible 
into an infected. Infected individuals spontaneously leave the infected state, reverting to susceptible. The well-mixed 
SIS system for sufficiently high infection probability possesses an endemic state, with a essentially constant level of 
infecteds, which however is subject to an exponentially small rate of extinction due to the existence of the absorbing 
state of zero infecteds. In addition, the deterministic rate equation is of Fisher-Kolmogorov type, and so a localized 
infection in the non-well-mixed case exhibits at the deterministic level an infection wave which propagates at constant 
velocity. However, the stochastic system exhibits not a bifurcation but rather a phase transition, characterized by the 
anomalous (for dimension D < 4) scaling exponents of the directed-percolation (DP) problem [7]. 

The connection between this complicated statistical behavior and the deterministic rate equations, which should 
be valid in the large- N limit (N being the total number of individuals, both susceptible and infected, on each site), is 
thus a natural topic for investigation, a study we initiated in a recent paper [5] (I). There we found that in one spatial 
dimension, the large- iV behavior was governed by a scaling law with an exponent which we called k ~ 0.66. For 
example, the phase transition point was shifted from the deterministic bifurcation point by an amount proportional 
to N~ K . Investigating the correlation length, £, we found that there was a scaling collapse so that £JV r , with r = 0.41 
was a function of N K times the distance to the deterministic bifurcation point. 

In this paper, we show that this behavior is in fact a nonequilibrium version of what one may call the Ginzburg 
crossover. A fundamental concept in equilibrium field theory is that of the Ginzburg criterion, which states under 
which circumstances the noise is relevant. This is of course what predicts the existence of an upper critical dimension 
(UCD), above which the noise does not affect the long distance behavior and so the scaling is mean-field like. The logic 
underlying the Ginzburg criterion implies that if one could "dial" down the noise, the system would look more and 
more classical, and a crossover (which we will call the Ginzburg crossover) between the classical and noise-dominated 
regimes should become apparent, with the noise dominating the very long distance behavior (below the UCD, of 
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course). This program has been implemented in the context of the equilibrium finite-range Ising model, where spins 
interact with all their neighbors out to a distance R. As R increases, each spin is interacting with what is more 
and more closely approximating the mean-field, and mean-field behavior at short scales sets in. The problem can 
be carried out analytically for the finite-range spherical model [5], and has been investigated via simulation in the 
finite-range Ising model [9]. 

Another fundamental epidemics model considered in I is the SIR infection model of Kermack and McKendrick [TU] , 
where recovered individuals are immune to further infection. The critical behavior of this model is governed by the 
dynamic percolation exponents, with an upper critical dimension of 6. In I we have carried a numerical investigation 
of the SIR model also, showing that the theory converges, in the large N limit, to its mean-field limit with scaling 
exponents k and r that differ from those of the SIS (DP) model. Here we consider again the SIR model and 
derive analytically the relevant exponents using the same theory of the Ginzburg crossover, applied to the different 
universality class of the SIR model. 

The main quantity used in the theory of epidemics to characterize the transmission potential of a disease is the basic 
reproduction rate Rq, which is the expected number of secondary cases produced by a primary case in a population 
that is completely susceptible [UJ. In the absence of demographic noise (e.g., in an infinite-dimensional model) the 
transition takes place at Rq = 1. Noise shifts the transition to higher values of i?o> but, as we will show below, the 
value Rq — 1 still admits a special feature: the renormalized distance between Rq = 1 and the actual transition point 
is N- independent as long as n < 1. This interesting feature allows one to examine the scaling properties numerically 
in a very efficient manner, as it saves the effort needed to identify the location of the transition point for each N 
separately. This feature is utilized here when we compare the expected results with numerical simulations. 



THE WELL MIXED SIS DYNAMICS AND THE TRANSITION ZONE 



We first review the well-mixed version of the SIS model, and the relation of the stochastic model to the deterministic 
equations that determine the evolution of the system. Although this stochastic model has already been analytically 
solved, the discussion allows us to present the concepts that we intend to use below and to set the mathematical 
framework used in the study of the spatial models. 

Let us consider a population of exactly N individuals, some of them are infected (I) and the rest are susceptible 
(S — N — I). The allowed processes are infection (with rate a/N, this is the type II model of [TT]) and recovery (with 
rate (3): 



S + I a -^2I I-t+S. (I) 

The corresponding master equation for the microscopic process can be formulated in terms of P n , the chance to have 
n infected individuals: 

P n = p {-nP n + (n + l)P„+i) + ^ [-n(N - n)P n + (n -l)(N-(n- l))P„_i] . (2) 
Defining (J) = ^ n nP n as the expected number of infecteds, one finds after index rearrangement: 

(/) = -m + ^J2 n ( N ~ n ) p - = (« - - ^< /2 >- ( 3 ) 

n 

The essence of the van Kampen f2 expansion is that this equation closes if (7) S> 1, so that the variance of I makes a 
negligible contribution, giving the standard logistic equation 

(i) = (a-/3)(I)-£L{I)*. (4) 

Since < (J) < N, it is necessary for N to be large, in order for the rate equation, Eq. Q, to be valid. This, 
however, is not sufficient. If a > j3, Eq. Q has an attractive fixed point at (I) = la = N(l — 1/Rq), where i?o = ot/fi 
is the primary reproductive number and Iq is indeed large if N is large, as required. Although the system admits 
an absorbing state at I = 0, the chance of a giant fluctuation that takes the system from Iq to zero is exponentially 
small in I , thus when N — > oo stochastic extinction (fadeout) is impossible once the system reaches its steady state. 
However, if the number of infected individuals in the initial state is small, stochastic effects are transiently present 
even in the N — > oo limit. For example, introducing one infected individual results in either short-time extinction 
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(with probability 1/Rq) or an endemic state (with probability (Rq — 1)/Rq). If Rq = 1 exactly, at the N — > oo limit 
the system performs an unbiased random walk in n, the number of infecteds, and the theory of first passage times 
tells us that the chance of extinction is still unity, but the probability P(q) to have q infection events scales like q~ 3 ^ 2 . 

At finite N the situation is more complex. Now the steady state of Eq. Q corresponds to a finite number of 
infected individuals in the endemic state, which mean that a finite, but large, fluctuation may cause a fadeout. The 
chance for such a fadeout is large when Rq is close to one, i.e., when the attractive fixed point corresponds to only 
a few individuals. Instead of having a sharp extinction to proliferation transition at Rq — 1, now the transition is 
"soft": defining A = Rq — 1 as the distance from the transition, Iq ~ NA; a metastable state exists only if this 
quantity (the distance of the stable solution from the absorbing state) is larger than the typical fluctuation size, y/N, 
thus a transition zone of width A ~ N^ 1 / 2 occurs between the extinction and the proliferation regimes. As shown 
in [T5], P(q) decays exponentially in the extinction phase A < 0, has a peak at exp(const ■ N) at the endemic phase 
A> 1/ VN, and decays like g~ 3 / 2 with a cutoff at N in the transition zone. Note that the width of the transition zone 
goes to zero as N approaches infinity, recovering the sharp transition at A = that characterizes the deterministic 
theory. 



THE ABSENCE OF SELF-INTERACTION 

The derivation of Eq. Q from Eq. ^ involves the neglect of 0(1/ N) terms. In particular one can easily see 
that the rate of infection when only one infected individual appears in a population of size N is a(l — 1/N), so the 
transition occurs at Rq = 1 + 1/N. This result reflects the most trivial effect of discretization, namely, the absence of 
self-interactions [T3]: an infected individual cannot infect itself, so the effective size of the population "seen" by the 
first infected is N — 1 instead of N. There are presumably other nonsingular 1/N corrections to the transition point, 
but for convenience we will refer to all these 1/N corrections as the "self-interaction" effect. 

Putting this fact together with the discussion of the last section, we realize that there are two N dependent functions 
that control the transition: one is the 0(1/ N) shift of the transition point, the other is the width of the "quantum" 
regime (the region above the transition point in which the system is controlled by demographic fluctuations) that 
scales, in the well-mixed limit, like N^ 1 / 2 . As N — > oo the shift is negligible with respect to the width of the transition 
zone, so there is only one scale in the problem, A ~ N~ x / 2 . However this behavior is not generic. As we will show 
below, in some cases the width of the transition zone is much narrower than 1/N, and in these cases one should take 
into account the two scales. 



SPATIAL SIS MODEL AND THE TRANSITION ZONE 



What happens if the system is extended? For the sake of concreteness let us focus on the example of an infinite 
one dimensional array of patches with N individuals on each patch. The probability per unit time of a susceptible 
on the nth site being infected by a given sick agent residing at this site is a(l — x)/N and of being infected by 
a given infected resident of one of the neighboring sites is a%/2A (in a d dimensional system, this chance will be 
ax/Nd). This corresponds to the "travelers model" considered in Ref. [13] ■ The deterministic rate, or mean-field 
(MF), equations, are 

L = -Pin + a ^~ X) I n (N - l n ) + ^(N- /„)(/„+! + I„-l) 

a X 2 , . a 2 ax 2 
= y" 1+ (a- /?)/„ - — /„ + — /V I (5) 

where V 2 stands for the discrete version of the Laplacian operator. The last, nonlinear diffusion, term, does not 
materially affect the dynamics (naive dimensional analysis shows that it is an irrelevant operator) . Without this term 
one recognizes, on the MF level, the celebrated Fisher (or FKPP [T5]) equation for invasion of a stable into an unstable 

phase, with a sharp transition at a = (or Rq = a//3 = 1), and front propagation with a velocity of 2y ax/3 A/2, 

since the effective diffusion constant is ax/2 and the net growth rate is a — f3 = /3A. 

What happen when stochasticity is taken into account? If N = 1, i.e., there is only one agent on any site and so all 
infections are nearest-neighbor (thus it is reasonable to take x — 1)j the stochastic process is known as the contact 
process, which undergoes a continuous phase transition from extinction to proliferation. The "effective" infection rate 
is smaller than a, since a sick agent cannot infect its neighbor if it is already sick. The transition happens at some 
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R c > 1, e.g., here for N = 1, R c « 3.297. While the exact value of R c is of course non-universal, the extinction 
transition, which belongs to the directed percolation equivalence class [7], admits three universal critical exponents: 

1. The spatial correlation length diverges as |A| _lyj -, where we introduce A = Rq — R c as the distance from the 
stochastic transition, as opposed to A, which measures the distance to the mean-field transition; in Id, v±_ ~ 1.09 

2. The temporal correlation length diverges like |A| _!y ii; in Id, vu s» 1.73 

3. Above the transition the steady state density of infecteds, Jo, grows like A' 3 ; in Id, j3 0.28 

The values of these critical exponents depend only on the dimensionality of the system and not on the microscopic 
details of the process. Above the critical dimension d = 4 the exponents take their MF values, v± = 1/2, un = 1. 
3=1- 

As N (the number of agents on a site) increases, demographic fluctuations become smaller. In the infinite N limit 
one recovers the MF transition described in Eq. First, the transition point moves back to R c — 1; second, the 
values of the critical exponent in this deterministic limit are equal to their MF values. For example it is clear from 
Eq. ([5| that above the transition the density scales linearly with A, i.e., that f3 — 1. Below the transition / is small 
and the nonlinear term in Eq. ([5| is negligible, hence if I(x, 0) = S(x), I(x, t) ~ exp{— x 2 /2Dt — At). The maximal 
density at x occurs when t ~ x/y AD; thus the spatial profile of total infections is proportional to exp(— x/£±) with 
£_L ~ 1/%/A, so that vf F = 1/2. 

At any finite N, though, close enough to R c the system is controlled by stochastic effects, as implied by universality. 
As N becomes large, the effects of stochasticity are restricted to a narrow region close to the transition point, which 
defines the width of the transition ("quantum") zone. 

In I, we have shown numerically that close to the transition point the spatial correlation length is given by: 

£l = AN- T (R C -R )-^ (6) 

where the transition takes place at R c = l + BN~ K . The values K ~ 0.66 and r ss 0.41 have been obtained numerically 
for different microscopic models that belongs to the DP equivalence class and seem to be identical for the different 
models up to the accuracy of the numerics. As long as k < 1, the region in the parameter space in which the system 
is controlled by stochasticity coincides with the interval between the stochastic and the deterministic critical points; 
i.e., it also scales like N~ K . Rescaling appropriately the correlation length and the distance from the transition, our 
numerics (see I) showed an whole scaling regime described by the function: 

n k-t/vx£-1/p± = f(AN K ) (7) 

The scaling function J- vanishes linearly at a positive value of its argument, which marks the transition point. Notice 
that what enters here is A = Rq — 1, so that the behavior at the classical transition point is controlled by the 
fluctuations, even though it is outside the range of the linear behavior of T . We will see later that the story is more 
complicated for k > 1. 

As discussed in the introduction, this scaling behavior is the result of a crossover between the deterministic theory 
and the critical theory as the critical point is neared. We will now use this to derive a scaling relation between r and 
k. Then we will obtain the value of k by calculating the Ginzburg criterion for the model. 

To connect r to k, one observes that the scaling function F{x) takes us from the stochastic regime at finite x (close 
to the transition) to the the deterministic regime at large negative x, corresponding to the region deeply below the 
transition. Even for \x\ large, the system may still be arbitrary close to the transition (A may be arbitrarily small) 
as long as N is large enough. This implies that in the x — > — oo limit, the correlation length must diverge like Ar 1 / 2 , 
independent of N. As a result the leading behavior of !F(x) at large negative x must obey F{x) ~ x 1 / 2u± . To cancel 
the N dependence in the expression 



one must have the scaling relation 



Given that we found k ~ 0.66, this implies a value of r w .40, consistent with our numerical findings. This scaling 
relation also implies that we can rewrite Eq. ^ as 

£ = iv K/2 t(an^ " x (10) 
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A similar argument is applicable to any of the quantities that diverge at the transition. One example that will 
be used below is the overall "mass" Mn of a cluster, namely the average total number of infection events before 
extinction. Utilizing the same scaling analysis, and the known mean field dependence Mm = 1/A, one expects that 
for re < 1, 



M N ~ N h 



g(AN K ^ 7 (ii) 



where Q vanishes linearly at the transition point, and <pM is the critical scaling exponent for the mass, 

7 « 1.24 (12) 



Eq. ( |11| ) is a useful relation that allows us to recover re directly from numerical simulations at fixed Rq. To demonstrate 
the critical exponents one has to locate exactly the transition point for any value of N; this is indeed a very tedious 
task. Instead, we can choose to simulate exactly at Ro — 1, which in our case implies Rq = 1. At this point the 
argument of the scaling function is exactly zero, independent of N, so the mass scales like N K . A plot of Mm/N k vs. 
N at Rq = 1 must converge to a constant in the large N limit. Below we will test this condition to verify numerically 
the predictions of our theory for re, as explained in the next section. 

However, this strategy works only for re < 1. As explained above, for higher values of re the trivial self-interaction 
shift of the transition point is not negligible in the N — > oo limit. Thus, as will be exemplified below, for dimensions 
where re > 1 one has to find first the transition point at R c = 1 + 0(1/N), and only near that point the transition 
region manifests itself. 



THE EXPONENT k AND THE GINZBURG CROSSOVER 

Determining re, thus, is enough to know everything about the quasi-deterministic regime. To find the value of re we 
adopt here a Ginzburg criterion approach, looking for the leading perturbative correction in inverse powers of N, and 
associate the stochastic regime with the region where this leading correction is 0(1). 

As a platform for the perturbative analysis we have chosen the Peliti-Doi field theoretic technique [T^] (see [T7] for 
details). Starting with the master equation for the SIS process, at a single site (zero dimensional system) with N 
individuals presented above. Eq. ([2| may be written as 

ijj = -Hi/) (13) 

where 

^ = ]TP». (14) 

n 

Using the creation-annihilation operators a\n) = n\n — 1) and a)\n) = \n + 1), the "Hamiltonian" takes the form 

H/P = (a 1 a - a) + i? (a f a - aVa) + -^(a f - l)(aWa). (15) 

Using the commutation relation [a^,a] = 1 and shifting from (that have a vacuum expectation value of unity |17j ) 
to a = — 1 one obtains: 

1 1 R 

H/P = (1 — -Ro(l — ^))^ a — ^o(l — ~~j^)o-Q- a ~ (aaaaa + 2aaaa + aaa) . (16) 

The first, "mass" term of the Hamiltonian determines the transition point: the system is in the active phase when 
the (renormalized) mass becomes negative. If N — > oo, an outbreak may occur at a > P, i.e., the transition happens 
when A = Ro — 1 = 0. The 1/N correction to this result reflects, again, the absence of self-interactions. 
Formally, the time evolution of ijj is given by: 

m = e- m ^{t = 0). (17) 

With the aid of time slicing and the coherent state representation one may arrive at a path integral representation of 
the evolution in time where the former creation-annihilation operators are replaced by complex-valued fields defined 
over a continuous space-time |17j : 

<tp(t)= [ VaVae- So( - a ' a) - Slia - a) iP(0) (18) 
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where 



with r7i = A 



S = / d d xdt a(x,t)[d t - £>V 2 - m]a(£,t). 



jf- and 



5 X = i? / d d x dt 



(1 )aaa (aaaaa 



2aaaa 



(19) 



(20) 



The renormalized values for all the constants in the problem may be obtained perturbatively by averaging over the 
cumulant expansion of exp{— S\) with weight exp(-So). The free propagator, in terms of spatial Fourier components, 



(a{k',t')a(-k,t)) = 8{k' + k)6(t - £') e (-fc 2 +™) (*-*•). 



(21) 



Here we are not really interested in the exact values of the perturbative corrections. All we are looking for is the 
width of the transition zone in the limit A — > and N —> oo. If a perturbative correction is proportional to N~ V A~ X , 
this correction becomes important (i.e., of order unity) when A = N~ y l x . There are many possible perturbative 
corrections with different x and y, but k is determined by the one that corresponds to the minimal value of y/x. In 
appendix 1 we will analyze the various elements of the perturbative expansion and conclude that 



d' 



(22) 



in particular k is 2/3 in one dimension, with almost perfect agreement with the numerical results reported in I. 
Moreover our result for a well mixed system (zero dimensions) is indeed n = 1/2, again with perfect agreement with 
the known results in that case. 

For 2d SIS, our expression predicts k = 1, so that the size of the stochastic regime is of the same order as the self 
interaction 1/iV corrections. As shown above, k determines also the relation between the average size of the epidemic 
and N when the infection rate takes its N — ¥ oo critical value, Rq = 1. Thus, in this case, we expect ~ TV. The 
data for this is presented in Fig. [T] The results are indeed consistent with the prediction; however the convergence is 
quite slow, much slower than in and 1 dimensions. 

In three dimensions, k > 1 and so the transition region is smaller than the 0(1/N) (self interaction) shift in the 
transition point. This leads to an interesting situation where there are two separate scaling regimes for large N. We 
will return to this point after first discussing the case of the SIR model. 



THE SUSCEPTIBLE-INFECTED-RECOVERED (SIR) MODEL ON SPATIAL DOMAINS. 

The other classic model of epidemics is the SIR model, which assumes that a recovered (R) individual cannot be 
infected again, so it is removed irreversibly from the "pool" of susceptible. The basic processes are: 

S + I ^ 21 J A it (23) 

The corresponding master equation for the microscopic process in a well-mixed population can be formulated in terms 
of to, the number of susceptibles, and n, the number of infected individuals: 

P„ = p {-nP m , n + (n + l)P m ,„+i) + ^ (-nmP m , n + (n - l)(m + l)P m+ i,„-i) . (24) 

In the deterministic limit, with the definition S — mP m n and I — ^ nP n ^ m and neglecting correlations (nm = 
nfh) one gets the equations: 

S = -~SI i=-pi+~SI k = pl, (25) 

where the last equation is just a consequence of the I dynamics. Since S = N — R — I , the two coupled equations, 
(again introducing A = a/ f3 — 1): 

i = f3AI-^I 2 -^IR R = fU, (26) 
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FIG. 1: The scaled "mass" of the aggregate, Mjv/iV, for the SIS model in two dimensions. The Ginzburg analysis suggests 
that, for large N, this ratio approaches a constant. Indeed, the plot shows that as N increased the ratio converges to a finite 
value. However this convergence is very slow, as implied by the iV~ 1/ ' 3 scaling of the abscissa. 



are enough to describe the system. The SIR dynamics does not support a nontrivial equilibrium steady state; instead 
at any site the epidemic disappears when t — > oo, leaving a finite density of recovered behind. This is manifested by 
the irreversible dynamics of R. 

Clearly, given I{x, t) one can solve for the number of recovered individuals at x: 



R(x,t) = /3 / I(x,r)dT, 



(27) 



Plugging that into Eq. (26 1 and adding terms that represent migration and discrete noise one gets: 

I = DW 2 I + (3AI - "-I 2 - °^rl [ I{t)<1t + jj(x, t)Vl 
N N J 



(28) 



where r\ is a delta-correlated noise, and D — a\/2 is the effective diffusion constant. Naive scaling analysis of Eq. 
(28) shows that the I 2 term is irrelevant and that the noise term becomes relevant when d < 6, as expected from the 
mapping to the dynamic percolation problem. Following [18| we integrate both sides of Eq. (281 from t = to oo, 
using J idt = and J °° I(t) J I(r)dtdT — V 2 [/o° I(t)dt] 2 , we arrive at 



DV 2 $ + - + ((x)V$ = 0. 



where &(x) = J °° I(x,t)dt. Note that the variance of the noise term in Eq. (29) must satisfy 



Noise 



dti 



i ^277(^)77^) v/y^/^y = / dt 1 i(t 1 ) = $ 

Jo 



(29) 



(30) 



justifying the form of the noise amplitude term in the <J> equation. Rescaling Eq. (29) by N we have (now m = /3A), 

1 



C(x)V$ 



(31) 



Eq. (31) may be analyzed perturbatively, as shown by [18] . by the same diagrammatic expansion used for the 
directed-percolation case (see Appendix 1) where the only difference is that the free propagator, instead of Eq. (21 ), 
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FIG. 2: The scaled "mass" of the aggregate, Mjv/TV, for the SIR model in D — 2, 3, and 4, showing the convergence to a finite 
value in the limit TV — > oo. 



is 



($(*)$(«/)>=<*(*' + *)_—. (32) 
fc z + m 

The first correction to the diffusion constant comes from the same self-energy diagram shown in Fig. [7J but here the 
correction is [TVA - ^-] -1 , thus n = 2/(6 — d). Accordingly, in both the SIR and SIS cases, we have that 

(33) 



d u — d 



where d u is the upper critical dimension. 

This result is consistent with the exact scaling of the transition region in SIR in dimensions, namely re = 1/3 [19j . 
It is also consistent with our numerical findings in I for the case of one dimension, where we found k « 0.41, to be 
compared with our prediction of 2/5. We can test our prediction for higher dimensions by again measuring the total 
mass at the classical transition point divided by TV K . This is presented in Fig. [2| The results are seen to converge 
relatively quickly to its finite TV = oo value in two and three dimensions, but show, similar to the SIS case in two 
dimensions, a very slow convergence in four dimensions. 



THE CASE OF n > 1 



As we have seen above, two scales are involved in the large TV limit. One is the shift of the transition point due 
to the absence of self- interactions, and this leads to 1/TV corrections for the critical reproductive number R c , and the 
other is the width of the "quantum" regime where fluctuations dominate the system behavior, the width of this region 
scales like N~ K . For d < d u — 2 we obtained n < 1 and the quantum regime is wider than the self-interaction shift, 
thus the effect of self- interactions is negligible. If d = d u — 2 both corrections scale like TV -1 and this leads to the 
slow convergence of the results to the large TV limit. We still have to consider the case where k > 1, i.e., where the 
quantum regime is narrower than the self interaction shift. 

For the SIS and SIR dynamics considered here, and for an integer number of dimensions, we have to consider the 
case k — 2 for d = d c — 1 (3d for SIS, 5d for SIR) and n = oo at the upper critical dimension. 

At d = d u the situation is trivial: k = oo means that the width of the transition zone is zero, since the system 
behaves (up to logarithmic corrections) like its mean- field (infinite dimensional) limit. Note the difference between a 



"Quantum" (fluctuation 
dominated) regime 



k< 1 



Ro=1 ^1+AN- 



"Quantum" regime 



k> 1 



Rq=1 Ro = 1 + AN- 1 R c = 1 + AN->+AN* 



FIG. 3: The two possible scenarios for large N scaling. If k < 1 (upper sketch) the 1/N self interaction shift is negligible 
with respect to the width of the quantum regime, thus the convergence to the deterministic limit is controlled by a single 
parameter N~ K . The case k > 1 (lower) is dominated by two scales: one that controls the distance of the transition point from 
its deterministic value, and the other that determines the width of the fluctuation dominated zone. 

well-mixed (Od) and the mean field (cod) cases: in the first there is a pronounced quantum regime at finite N. In the 
second each point has infinite number of neighbors so the "effective N" is infinite even if the number of individuals 
at each point is finite. 

What remains is d = d u — 1, namely three dimensions for SIS and five dimensions for SIR. In these cases, n = 2, 
so the "quantum" regime has a very small width (of order N~ 2 ) around the quantum transition point, which in turn 
is at a much larger distance (of order 1/N) away from the deterministic transition point Rq — 1. The situation is 
summarized in Fig. [3] 

Although the transition point converges to Rq = 1 in the infinite N limit, this convergence is slower than the rate 
in which the quantum zone shrinks around this point. This gives rise to two different scaling regimes, one of width 
1/N and the second of width 1/N 2 . We can see this behavior, again, by studying Mm, the total mass of the infection, 
now as a function of Rq, the "bare" reproductive number. 

In the outer region, with width 0(1/N), the mass obeys the scaling law 

M N = ATm Gout (Aa) . (34) 

Now, far from the transition point (say for fixed Ro slightly below 1), at large enough N the dependence of Mm on 
this distance must approach its MF limit, Mm ~ 1/(1 — Ro), independent of N. This implies that for large negative 
argument, G(x) ~ —1/x, and that r^ 1 ' = 1. Since the transition point is at R c « 1 + (A/N) + (B/N 2 ), where A and 
B are some constants, Mm must get large as AN approaches A. Since in this outer region, fluctuations are small, 
the incipient divergence of M is mean-field like, so Gout diverges as Gout ~ c/(A — x), so that for Ro near, but not too 
near R c , Mm behaves as 

C 

M N « n (35) 

1 + A/N-Ro V ; 

This behavior is demonstrated in Fig. |ij where N/Mm is plotted versus A7V, for SIS in three dimensions in the 
upper panel and for SIR in five dimensions in the lower. We see that there is a very slow convergence to an asymptotic 
curve. This slow convergence to the asymptotic scaling limit is reminiscent to what we encountered in the case of 
d = d u — 2. The large N line is straight, but does not converge to zero at Ro — 1, since the actual transition happens 
at Rq ~ 1 + A/N . Although at large N, the distance of Ro = 1 from the transition shrinks to zero one observes no 
"quantum" effects in the outer region since the width of the quantum regime shrinks even faster. 

As we approach very close, of order a small fraction of 1/N 2 , to the phase transition point, the fluctuations become 
significant and Mm diverges as 



M N =AN- T ™{R c -Ro)-^ 



(36) 
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FIG. 4: Upper Panel: The inverse of the scaled "mass" of the aggregate, N/Mn, for the SIS model in d = 3 as a function of 
N(Rq — 1), for various N. The behavior for large N(Rq — 1) is consistent with Mm = 1/(1 — Ro)- The data labelled N = oo was 
obtained by fitting a quadratic curve in iV~ 1//3 to Mn for fixed N(Rq — 1) and extrapolating. This N = oo curve fits well to 
Mjv = N/(7.2 — N(Ro — 1)), corresponding to a shift in the critical Ro by an amount 7.2/ TV '. Lower Panel: The inverse of scaled 
"mass" of the aggregate, N/Mjsr, for the SIR model in D = 5 as a function of N(Rq — 1), for various N. The behavior for large 
N(Rq — 1) is consistent with Mn = 1/(1 — Ro)- The data labelled N = oo was obtained by fitting a quadratic curve in TV -1 / 3 
to Mjv for fixed N(Ro — 1) and extrapolating. This N — oo curve fits well to Mn = N/ (47.6 — N(Ro — 1)), corresponding to a 
shift in the critical Ro by an amount 47.6/iV. 
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where 7 is the scaling exponent for the mass, which for DP in three dimensions is 7 ~ 1.24 and is approximately 
1.2 for percolation in five dimensions. The general scaling law in the inner region, of width 0(1/N 2 ) is then 

M N = N T % [g m (AN 2 )]~\ (37) 

where Qi n {x) vanishes linearly at x = 0. For large negative argument, this has to match onto the outer behavior for 
AiV <C 1. This is possible if Qi n {x) ~ — Cx as x — > 00 and r^" = 2. 

Accordingly, the plot of (iV _2 Mjv) -1 / 7 vs. AN 2 shows the inner scaling function Qi n in the large N limit and goes 
linearly to zero at the transition point. This behavior is demonstrated in Fig. [5] for both the 3d SIS (upper panel) 
and the 5d SIR (lower panel) models. 



FRONT VELOCITY 



In the wake of Brunet and Derrida's [3] pathbreaking work on the large N behavior of the front velocity in Fisher- 
type systems, there has been an enormous amount of attention devoted to this issue, including a rigorous proof of the 
original heuristic arguments. It is thus natural to ask how this work relates to our current findings. The first thing 
to note is that the limits addressed here and the result of [3] are different. The Brunet-Derrida limit corresponds in 
our language to fixed A, N — > 00, whereas we are interesting in the limit A <C 1, N ^> 1, AN K ~ 0(1). 

We first investigate the behavior in the immediate vicinity of the transition point, restricting our attention to the 
Id SIS model. In the immediate vicinity of the transition point, both the spatial correlation length, and the time 
correlation scale, £y diverge. It is expected then that the velocity will scale, in this regime, as the ratio of £j_ to £y: 



B_ L N~ n ^- L ~^A~ v - L 
5iiiV~ K(,y ii~ 1) A~ I/ ii 



DA"!!-^ = B V N~ K/2 (AiV K p 



(38) 



Since > z/j_, the velocity vanishes as the transition point is neared, just as in the classical theory. Furthermore, 
V\\ — u± — 1/2 > 0, so the velocity increases with N for fixed A. This is consistent with the Brunet-Derrida asymptotic 
result, which also has the velocity rising with TV at fixed A. 

In the classical limit the front velocity is given by v ~ V^A, independent of N. One is tempted, then, to write, 
in analogy with our other scaling laws, v s» N~ K / 2 'H(AN K ). The problem with this is that, while in the continuum 

classical limit, the velocity is proportional to vA, on the lattice this is true only for small A. To work with discrete 
agents and to define their local density one should implement some UV cutoff, so even for off-lattice models the 
relevant result is the one obtained for a lattice. The classical lattice velocity vl satisfies the equation (see |22| ) 




Svl\ = £A 
\axj ax 



(39) 



so that, for large A, the velocity grows as A/ln(A), as opposed to 
relation for v, it is preferable to find a scaling relation for 



7(«) 



A. Thus, instead of trying to find a scaling 




1/2 



(40) 



which, for v — vl, is precisely equal to V A. In Fig. |6j we show the scaling collapse of g(v)N K / 2 versus AN K . The 
Brunet-Derrida effect, namely the anomalously slow approach to the classical velocity, is apparent from this graph, 
where even for AN~ K ~ 60, the scaling curve is very far below the classical result. 
In more detail, for large positive argument, the Brunet-Derrida result implies that 



H(x) 



1 



9ir 2 



4 In 2 a: 



(41) 



This corrected classical result is also show in Fig. |6j where we see quite good agreement, especially considering the 
relatively small values of N involved, compared to those necessary to achieve even semi-quantitative agreement with 
the Brunet-Derrida correction at A ~ 0(1). 
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FIG. 5: (color online) Top: The scaled "mass" of the aggregate, (M^/N 2 ) 1 / 7 , as a function of the scaled "inner" variable 
N 2 (R C — Ro), for the SIS model in D = 3. We used the value 7 = 1.236. Bottom: The scaled "mass" of the aggregate, 
(Mn/N 2 )~ 1 /~' , as a function of the scaled "inner" variable N 2 (R C — Ro), for the SIR model in D = 5. We used the value 
7 = 1.23. 
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FIG. 6: The scaling collapse of the scaled transformed velocity, g(v)N K ^ 2 , where g(v) is given in Eq. (40 1, versus AiV" for the 
Id SIS model, with /3 = 1, x = 0.2. The "classical" result is (R — l)iV K , whereas the "corrected classical" result is given by 
Eq. pll. 



SUMMARY AND DISCUSSION 



Along this paper we have studied, numerically and analytically, the convergence of the stochastic process to the 
deterministic rate equations when the number of particles is large. In spatially extended model there are two param- 
eters that control the convergence: the number of particles per site N and the distance from the transition point. 
Together, these parameters yield a region of size AN~ K above the phase transition point; within this region the system 
is dominated by demographic noise and the deterministic equations fail to describe it accurately. 

The value of K has been found before using an extensive analysis of zero dimensional [T2J []j|J US H3] and one 
dimensional [5] models. It turns out that this particular result may be derived directly, for any dimensionality using 
the Ginzburg analysis. For the fundamental models considered here it depends only on the difference between d u , the 
upper critical dimension, and d, via k = 2/{d u — d). 

Clearly, this general expression stems from the fact that the leading perturbative correction (i.e., the diagrams that 
lead to an infrared divergence in the highest dimension, which is thus the upper critical dimension) is proportional 
to 1/iV, since it involves an average over two noise terms, each is proportional to 1/y/N. This seems to be a generic 
property of stochastic processes and will be interesting to find out a model for which this general argument is not 
applicable. 

Below d u — 2 k < 1, and the self-interaction shift is negligible at large N. In this case the point Rq = 1 is 
peculiar: its normalized distance from the critical point (the distance divided by the width of the quantum regime) 
is N independent. Accordingly, the divergence of various observables at this point is determined solely by N K . This 
feature facilitates the numerics, since one can extract the value of the exponent without finding R c . If k > 1 this is 
no longer the case, and to locate the quantum regime one has to first identify the transition point. 

Although the SIS and SIR processes serves us here as an archetypic stochastic processes that belong to the most 
pronounced equivalence classes of out-of-equilibrium transitions, they are also interesting models for epidemiologists. 
Several attempts have been made, recently, in order to understand better the role of fluctuations in individual- 
based, spatially structured epidemic models. The results presented here practically solve this problem for the case of 
subpopulations on a lattice considered in [24]. 

In the common case of zoonotic infections the pathogens first emerged from animal reservoirs, inducing a "stuttering 
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FIG. 7: Element of the diagrammatic perturbative expansion. The terms that appear in Eq. (20 1 (upper line), the self energy 
diagram £ and the 1-loop correction to the three point vertex (lower part). 



transmission" stage in which R < R c , and reaching the phase of sustained transmission (human outbreak) only due 
to pathogen evolution (in human environment) to Rq > R c |25j . If Rq is growing slowly to larger value (as opposed 
to a major evolutionary step caused by a single mutation) the pathogen must cross the quantum region, where the 
size of the outbreak (the number of infections, and hance the chance for the next evolutionary step to occur) is 
simply Mn(R). With an appropriate knowledge about the adaptation process of the pathogen, it will be quite easy 
to implement our results to obtain the chance of an outbreak. 
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APPENDIX 1 



Here we show some elements of the perturbative expansion of the action (19 20 1 and the terms that determine the 
leading correction for large N, as explained in the text. 

The elementary diagrams that appear in the perturbative expansion are shown in the upper part of Fig. [7] Of 
those, the first two appear in the Reggeon field theory and yield the one- loop renormalization of the mass and i?o- 
The diagrams involved are presented in the lower part of Fig. [7J 

With the bare propagator, Eq. (21 1, one can see that the leading correction to the mass behaves like 



1 

N 



dq 



This implies that q scales like y/m and hence close to the transition the result is proportional to A~( 2_d )/ 2 , thus from 
this diagram one would get k = 2/(2 — d) (x = (2 — d)/2, y = 1, see text). The triangular diagram that provides the 
correction to the coupling constant scales like 



1 

iV 2 



dq 



(q 2 



so it corresponds to k — 4/(4— d). However, the corrections to the diffusion constant are given by the second derivative 
of the self-energy diagram with respect to the incoming momentum, and this contribution is proportional to 

1 f q d - x dq 
N J (q 2 + m) 2 ' 
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and this term yields the minimum value k = 2/(4 — d) given in Eq. (22 1 
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